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ON THE CONSERVATION AND CONVERGENCE TO WEAK SOLUTIONS OF 

GLOBAL SCHEMES* 

MARK H. CARPENTER*, DAVID GOTTLIEB*, AND CRT- WANG SHU§ 


Abstract. In this paper we discuss the issue of conservation and convergence to weak solutions of several 
global schemes, including the commonly used compact schemes and spectral collocation schemes, for solving 
hyperbolic conservation laws. It is shown that such schemes, if convergent boundedly almost everywhere, 
will converge to weak solutions. The results are extensions of the classical Lax-Wendroff theorem concerning 
conservative schemes. 

Key words, conservation laws, conservation, weak solutions, convergence 
Subject classification. Applied and Numerical Mathematics 

1. Introduction. We are interested in numerical solutions to the conservation laws: 


(1.1) u t + f(u) x = 0, u(x,0) = u°(x), -1 < x < 1. 

Here we have written (1.1) in the one dimensional form, but the results of this paper are also valid for multi 
dimensions. 

The purpose of this paper is not to study the issue of convergence. We actually assume that the numerical 
solution converges boundedly a.e. (almost everywhere), to a certain function u(x,£). More precisely, for a 
numerical scheme defined at the (uniform or nonuniform) grid points Xj,0 < j < TV, with Ax = max(x J+ i - 
xj) and Vj (t) as the numerical solution at x — Xj, we define the function uaj-(£, 0 by 

(1.2) WA x{x,t) = Vj(t, ), Xj < X < Xj + i, 


and assume that ua*(x,£) is uniformly bounded with respect to x, and Ax, and, as Ax — » 0, v& x (X)t) 
converges pointwise a.e. to u(x,£). See, e.g., [5, 16, 17, 20] for discussions, in the scalar case, of convergence 
of some of the schemes studied in this paper, under the L°° boundedness assumption. We will concentrate 
on the issue of whether the limit function u(x, t) is a weak solution to (1.1), that is whether it satisfies 


(1.3) 


a: 


{u{x,t)(j>t{x,f) + f{u(x y t))(j). 


■ dxdt- J u°(x)<j)(x,0) dx = 0 


for any smooth function <£(x, t) which is compactly supported. Also, in this paper we only consider semidis- 
crete method-of-lines schemes, i.e. schemes which are discretized in the spatial variable(s) only. 

The classical result in this area is the famous Lax-Wendroff Theorem [11]: 


Theorem 1.1. (Lax and Wendroff) If the numerical solution of a conservative scheme: 

C- 4 ) + 
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where the numerical flux 

(1*5) /j-fl ~ / ( v j~p> •••) v j+q) 

is local (i.e. p and q are constants independent of Ax), Lipschitz continuous in every argument, and consistent 
with the physical flux /(u, — f(v), converges boundedly a.e. (almost everywhere) to a function u(x> t), 

then u(x,t) is a weak solution to (1.1). 

■ 

The proof follows easily from a summation by parts and an application of the dominant convergence 
theorem. See LeVeque [13] for a slightly different version of this theorem and its proof. 

The Lax-Wendroff Theorem, however, does not cover global schemes, i.e. schemes which can not be 
written in the form (1.4) with a local flux /j+i- Examples of global schemes include the compact schemes 
[12, 2, 3, 5], and various spectral Galerkin or collocation schemes (Fourier, Legendre, Chebyshev) [9, 1, G]. 
We will extend the Lax Wendroff Theorem to these global schemes in this paper. 

We remark that there are discussions in the literature about schemes which are not of the conservative 
form (1.1) but nevertheless still converge to weak solutions. One such example is the class of schemes for 
general curvilinear coordinates, see [19] for a proof that such schemes actually do converge to weak solutions. 

In [8] the authors discussed conservation issues of Chebyshev methods. However, they only considered 
the mean of the solution, that is, they verified that the limit solution satisfies (1.3) with <f>{x,t) = 1. 

The organization of the paper is as follows. In Section 2 we discuss compact schemes. Section 3 contains 
the Legendre collocation method, while Section 4 discusses the Chebyshev method. Section 5 discusses the 
Legendre approximation in the multi-domain case. 

To close this section we mention that in this paper C (or c) is a generic constant. 

2. Compact Schemes. Compact schemes are methods where the derivatives are approximated by 
rational function operators on the discrete solutions. We consider compact schemes defined on a uniform 
grid XjsO < j < N. For example, a fourth order central compact approximation to the derivative is [12]: 

(2T) - ((v x )j-i + 4 (v x )j + (v x )j+i) = - (u/+i — vj_ i) 

and a third order upwind compact approximation to the derivative is [5]: 

( 2 * 2 ) ^ {-(Vx)j - 1 + ${ v x)j - (v*)j+ 1 ) = “ (3 vj - 4uj_i 4- vj - 2 ) . 

Adequate boundary conditions must be used for the compact schemes, to retain accuracy and stability, see 
[2], [3] for details. Together with boundary conditions, a compact scheme for (1.1) can be written as 

(2.3) Pvt + QJ{v) — -t(v b — 9b) 

where v = (u 0j is the numerical solution, r is a constant, v B = (u 0 , 0, ..., 0 } v N )* is the boundary part 

of the numerical solution, and gs = (<?o, 0, •••, 0, 9n)* is the given boundary data. Depending on the wind 
direction, one or both of the first and last components of vb and 93 may also be zero(s). The matrices P 
and Q satisfy the following conditions [3]: 

• P is symmetric, and satisfies 

(2.4) P<j)-<j)~ O(Az) 

Here and below, <j> = (<f>(x 0 ), ..., <j>(x at))* , and <j>{x) is an arbitrary smooth (C 1 or smoother) function. 
O(Ax) for a vector means that each component is bounded by a constant times Ax, and the constant 
depends only on the derivatives of <t>{x)\ 
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• Q is “almost” anti-symmetric, that is: 

(2.5) Q + Q* — R + S 

where R = (r*j), and = 0 except for r 0 o and 5 is either identically 0 for the central compact 

schemes, or satisfies 

(2.6) Set) = O(Az) 

for the upwind compact schemes, where <p is defined as before. Also, Q is at least a first order 
approximation to the derivative: 

t (2.7) Q<S>-4>x =0(&x). 

We can easily verify that all the compact schemes in [12, 2, 3, 5] satisfy the above conditions for P and 

For such compact schemes we can state the following proposition: 


Proposition 2.1. If the solution of the compact scheme (2.3) converges almost everywhere to a function 
u(z,£), then u(x, t) is a weak solution to (1.1). 

Proof: For any compactly supported, C 2 function <£(;r,f), we denote by <j> = (<£(# 0 , t), <f>(%, £))*> ^ft 
multiply (2.3) by cj )* , and integrate over [O.T] to obtain: 


f 


<j ) * {Pvt + Qf( v)) dt = 0 


due to the zero boundary conditions of 0. Now integrating by parts in t for the first term, taking a transpose 
of the equation (which is a scalar), and using the symmetry of P and condition (2.5) of Q, we obtain: 

rp rp 

- [ (v'Pfa + f{v)*Q<t>) dt - {v*P<p) | t=0 = - f f(v)*s<f>dt 

Jo Jo 

Or, considering (2.4), (2.6), (2.7), and the uniform boundedness (with respect to the mesh size A:r) of v , 

(2.8) - f T (v*4H + dt - Uo = 0(1) 

Jo 


where the constant term 0(1) results from a summation of N = 0 ^ r) terms of O(Ax) quantities; the 
constant depends only on the derivatives of <j>{x ). 

Recalling the definition of the function v^ x {x^t) in (1*2), we can multiply (2.8) by A.r to obtain 

(2.9)-/ J (v& x (x,t)<j) Ax {x,t) + f{vAx{x,t))<j> Ax {x,t)) dxdt - J v Ax (x,0)<j> Ax {x,0) dx = 0{Ax) 


where 


(2.10) <j> Az (x,t) = <f>{xj,t), <j> Ax {x,t) = 4> x (xj,t), (p Az (x,t) = <p,(xj,t), xj<x<x j+ 1 , 

By assumption, v& x (x, t) converges to u(x,t) boundedly a.e. There is no problem about the uniform 
convergence of <t>^ x (x,t) and <^ x (:r,t) due to the smoothness of 0. By the dominant convergence 

theorem, taking the limit as A# — > 0 in (2.9), we obtain (1.3). This proves that u(x , t) is a w r eak solution of 

( 1 . 1 ). 
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3. Legendre Spectral Collocation Schemes. The Legendre collocation method can be written in 
the following way: 


(3.1) 


du N {x,t) dI 2 \f(u N (x,t)) 
Ot dx 


SV(u N {x , <)) + Bu N {x,i) 


w r here w/v(x, £) is the numerical solution which is a polynomial of degree at most N in x, I ^ is the Legendre 
interpolation operator, i.e. for any function g(x ), Isg{x) is the unique polynomial of degree at most N 
satisfying Ij\rg(xj) = 9( x j) the N + 1 Legendre Gauss-Lobatto points Xj y which are the zeros of the 
polynomial (1 — x 2 )P Ni where P/v is the Legendre polynomial of degree N. 

The term SV is the spectral viscosity term needed to stabilize the scheme and in order for the assumption 
u v& x {x^t) converges boundedly a.e. to a function u{x , £)” to be realistic. We consider here the superviscosity 
term 


(3.2) sv(<w) = 

6 is the superviscosity coefficient, s is an integer growing with N [21, 10, 14, 15]. We remark that this 
superviscosity term is equivalent in practice to a low* pass filter. 

Finally, the boundary term Bujy{x 9 t) could be either 0, or 


£<* 


* 2 ) 


d_ 

dx 


u N (x, t ) 


t(ujv(M) -ffl)(l + x)P' N (x), 


or 


t(u n {- 1,0 — <7- 1 ) ( 1 -x)P' N (x), 

or a combination, depending on the wind directions at the boundary points. Here r is a constant chosen for 
stability and g\ and g-\ are functions of the time only. 

Let <j>{x,t) be a test function in Cq° . Tak then clearly are polynomials of 
degree at most A T -1 and vanish at both boundary points x — dbl. Also <pN-\ ( x , t) -> t) y ((j>M-i)x{Xi t) 
4> x (x, t), and -> <f> t (x,t) uniformly. 

We denote now’ by 

C/i 9) = J f(x)g{x)dx 

and by 

AT 

(/>P)n = 53 f(x j )g{x j )uj j 

k=0 

where u>j > 0 are the weights in the Gauss-Lobatto formula. We note that ( f,g)n = (f,g) if fg is a 
polynomial of degree at most 2A r — 1. 

We first show that the boundary terms do not cause a problem: 

Lemma 3.1. 

(3.3) (6 n _ u Bu n ) = 0. 
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Proof: We start by observing that 


(4> N - U Bun) = {4>n-i,Bu n ) n . 


Bun vanishes for the inner Gauss-Lobatto points and (f>N-i vanishes at the boundaries and therefore the 
Lemma is proven. 

■ 

With (3.3), we multiply (3.1) by <f>jv- i(x, £), integrate over x , and integrate by parts for the second term 
to obtain 


(3.4) 


1 d<j>u-\{x,t) 




d_ 

dx K 


Ox 


un(x, t) dx. 


We now estimate the right hand side of (3.4): 

Lemma 3.2. 


(3.5) 


e(— 1 )* rl 
N^‘oo N 28 ~ 


lim 




4>n-\ (x, t) 


Ox dx 


un{x , t) dx = 0. 


Also, the quantity under the limit sign is uniformly bounded with respect to t . 
Proof: Since <j> n-i is a polynomial of degree N — 1, 

v-i 


<t>N-i(x,t) = <f>k,N(t)Pk(x) 


k-0 


where (j>k y N{t) are the collocation Legendre coefficients of the test function (ft. Note that 

p k (x) = (-iyk’(k+ \yp k (x). 


-a-T 2 )— 

dx ( dx 


and therefore 


(4 > n -uSV(u n )) = (-iyj^Y^k%k + iyt k>N (t)u k>N (t)(p k ,p k ), 


k = 0 


Here Uk,N are the Legendre collocation coefficients of uyv- We note that as a consequence of the uniform 
boundedness of un(xj) and the fact that 4> is in Cg°, 

(3.6) < p** 

This implies (3.5) and the uniform boundedness of the quantity under the limit sign with respect to t. 

m 

We thus only have to deal with the left hand side of (3.4). We integrate (3.4) in t, integrate by parts for 
the first term, and use Lemma 3.2 to obtain: 


L L 


t) 


0<f>N-i (x,t) 


dt 


+ I N f{u N (x,t)) 


d^N- iQM) 

dx 


) 


dxdt 


(3.7) 


- J u N (x,0)<}> N -i(x,0)dx = o(l). 
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It looks like we can immediately take the limit as in the Lax-Wendroff Theorem. The trouble is that, we 
have only assumed the uniform boundedness of hence of but this does not imply the 

uniform boundedness of either u N (x, t) or I N f(u N (x, t)) due to tlie lack of regularity. 

We need the following Lemma: 

L emma 3.3. Let be the piecewise linear polynomial taking the values uw(xj,t), then 
(3 8) (u N (x, t), = ^ + o(1); 

0-9) (l N f(u N (x,t)), = (/n/(^,), ^~‘ (M) ) +o(l). 


Proof: We will switch back and forth between integrals and quadrature summations: 


■where 


u N (x,t), 


dt 


^ ( t o 0<pN-l(x,t)\ 

j = 

N 

E , .■.d<pN- 1 (x.j,t) 

U N {xjJ) U)j 

J = 0 Ul 

= J 1 6(x)vat {x, )f i (i, t)dx 


0(X) = 


Xj+1 - Xj 


Xj ^ X ^ Xj j . 


In [18] it has been established that 9{x) is uniformly bounded and converges a.e to 1 as Ax — > 0. 
This proves (3.8). The proof for (3.9) is similar. 


We can state now 


Theorem 3.4. If the function v± x (x,t) defined in (1.2), obtained from the solution of the Legendre collo- 
cation scheme (3.1) at the Legendre Gauss-Lobatto points xj, converges almost everywhere to a function 
w(x,t), then u(x,t) is a weak solution to (LI). 

Proof: By assumption, v& x (x,t) converges to u(x, t) boundedly a.e. Also, there is no problem about the 
uniform convergence of t) and (0A r -i)^ x (£, 0 due to the smoothness of <p. Using 

the dominant convergence theorem, taking the limit as Ax -» 0 in (3.8), (3.9) and (3.7), we obtain (1.3). 
This proves that u(x,t) is a weak solution of (1.1). 


We close this section by commenting on other spectral viscosity terms in (3.1) that stabilize the Legendre 
method. One such term is 


ox 


du N 

dx 


where the spectral viscosity operator Q is defined by 


N 

Q4> - ^Q k 4>kPk{x) 

k = 0 


G 


where 


N 

<f, = Y,kPk(x) 

A=0 

and 


1 > Qk > l - 


Qk = o ,k < mu 


with mjv growing witli N . 

We can establish also for this viscosity term that 


(0AT-1,SI (uat)) -> 0 


and therefore the result above holds also for this kind of spectral viscosity. 


4. Chebyshev Spectral Collocation Schemes. In this section we consider the Chebyshev collocation 
schemes. These are more difficult to analyze than the Legendre method because of the weight function y= ==^ . 
The Chebyshev collocation method can be written in the following way: 


(4.1) 


du N (x,t) dJ A f (u N {x,t.)) 


dt 


dx 


e(-l) 9 

IV 2 ®" 1 


ox 


2s 


u N (x,t) + Bux(x,t) 


where again u A '(x,t) is the numerical solution which is a polynomial of degree at most N in x, Jjv is the 
Chebyshev interpolation operator, i.e. for any function g(x ), Jjvg{x) is the unique polynomial of degree at 
most N satisfying Jn9{xj) = g(xj) at the A r + 1 Chebyshev Gauss-Lobatto points xj. e is the superviscosity 
coefficient, s is an integer growing with N [21, 14, 15]. We remark again that this superviscosity term, which 
in practice is equivalent to a low pass filter, or a similar vanishing viscosity term [1G, 17], is needed in order 
for the assumption u v& x (x,t) converges boundedly a.e. to a function u(x,t) v to be realistic. Finally, the 
boundary term Buw(x,t) could be either 0, or 


r{u N (l, t)-gi){l + x)T' n (x), 


or 


r{u N {-l,t) - 5-i)(l - x)T' n (x), 

or a combination, depending on the wind directions at the boundary points. Here r is a constant chosen for 
stability and g\ and g-i are functions of the time only. 

Let <j>(: r, t) be a test function in Cq , that is, all x derivatives of </>(x, t) up to order 5 vanish at the boundary 
points x = ±1. Such test functions are, of course, dense in Cq. It follows that (1 - x 2 )~ 3 / 2 <fr(x, t) is in Cq- 
We denote the (N — 5)-th degree Chebyshev interpolation polynomial of the function (1 - x 2 )- 3 / 2 <fi(x, t) by 

(4.2) £/v-5(z, t) - Jtf_5((l - x 2 )~ 3/2 4>{x, t)), 

and note that 

Cv-5(x,f) -»• (1 - X 2 )~ 3/2 <t>{x,t), 

dtN-6 (X,t) 0 ((1 - x 2 )~ 3 / 2 <j>(x,t)) dtK-f>(x,t) d{(l-x 2 r 3 ' 2 <i>(x,t)) 

dx dx ’dt dt 
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uniformly. We now take 
(4.3) 


Vbv-i (M) = (1 -X 2 )-Z N -a{x,t), 


then V'jv-i is a polynomial of degree at most N — 1 and vanishes at both boundary points x = ±1 together 
with its first, and second x derivatives. Moreover, it can be easily verified that 


(4.4) 


$n- iOM) 
yi — X 2 




v/T 


!?).—<*•'»■ (W), 


(x, t) , 


uniformly. 

We again first show that the boundary term does not cause a problem: 

Lemma 4.1. 

(4-5) j 1 ^+x)r N (x) '^t^- dx = 0, j\l-x)r N (x)^- l ^dx = 0. 

Proof: We only prove the first equality. Zero boundary values of V^v- i(x, t) and its first x derivative imply 

/_', (1 + *")W * “ - r »<*> (^7r^" + *>). * 


=/: 


= 0. 


Tn(x) f ctyjy_i(z,£) 

VT^Z 2 V 0 * 


(1 + x) + 


^.v-](j-.Q 

1 - X 


dx 


The last equality is due to the fact that 


dll>N-l(x,t) „ rp N _ l(x,t) dll>N-l{x,t) „ , x , ,, , X/1 2x ,. , JN 

-(1+*) + — 5 — = HI (1 + z) + (l+x)(l - x-)t N _ 5 (xj) 


dx 


1 - X 


dx 


is a polynomial of degree at most A 7 - 1, hence is orthogonal to T N (x) with the weight 

■ 

With (4.5), we can now multiply (4.1) by , integrate over x, and integrate by parts for the 

second term to obtain 


/ 


1 lJ>N-l(x,t) du N (x,t) 


Qt 


(4.6) 


dX ~ f t ( ^1 - J Nf{u N {x,t))dx 

€(-!)* f l *l>N-l(x,t ) 

n 2s ~ i J_y 


(l-x 2 ) 


Ujv(x T t) dx. 


We now estimate the right hand side of (4.6): 

Lemma 4.2. 

e(-l) s f l i!> N -i(x,t.) 


(4.7) 


liin 


jr r 
s - 1 7-i 


\/l - x 2 -^- 

ox 


2s 


uat(x, £) dx = 0. 


AT->0C A r2s_1 v /fZ“^2 

Also, the quantity under the limit sign is uniformly bounded with respect to t. 

Proof: Integrating by parts 2 s times, and noticing that the boundary terms are always 0 because of the fact 
that ^tv~i(x, f) vanishes at the boundaries and because of the factor y/1 — x 2 , we obtain 


<«> cm 


t) 


v/l - 

ox 


u N (x 


:\ t) dx = J 


1 uyy(M) 

V Y— x 1 


yj\ — x : 


d_ 

dx 


-| 2 8 


^v_i(x, t) dx. 


8 



Recalling the definition of £/v_ t5 (x,f) in (4.2), we have 

A T -5 

& v - 5(*,0 = E &(<) r fc (*) 

A=0 


where £*.(/) are the collocation Chebyshev coefficients of the Cq function (1 — x 1 ) 3 / 2 <^>(x,f), hence 

(4-9) 16(01 < ' 

Now, by the relationship between t/W-i and £a'-s in (4.3): 

A T — 5 

1 p N -l(x,t) = (1 - x 2 ) 2 E ik{t)T k {x) 

k=Q 

1 TV— 5 

= 3 (1-T 2 (x)) 2 E6(0 ^(x) 

fc =0 

A T — 5 A T — 5 

= i E www - 5 E & w ( t h2W + T k - 2 {x)) 

k-0 k=0 

1 A r — 5 

+Tg E (^+ 4 ^) + mix) + r*_ 4 (x)) 

° A— 0 

- A r — 1 

= 7h X (0 + 6£fr(t) - 4^fc + 2(0 + £fc+4(*)) 7*0*0 

° fc =0 
A r — 1 

= XI ^*w T *o*o 

*=0 


where we take the convention that £*(£) = 0 for k < 0 or k > N — 5. This, together with (4.9), clearly 
implies 

(4.io) hMOl < 

We now 7 use the equality 


(4-11) 


y/l — x 2 


d_ 

dx 


2s 


^JV-lOM) 


N — \ 


E(-l ) s k 2 °Mt)T k (x) 

k = 0 


and the integral-quadrature equivalence: 


/ 


-i Vi — a -2 L 




2s 


^AT_i(x, t) dx 


AT 

j— o 



>A 

ax 


1 2s 


^_i(x,0 



where xj and wj are the nodes and weights of the Chebyshev Gauss-Lobatto quadrature formula, because 
the integrand 


U N (x,t) 




2s 


tl> N -i(x,t) 


is a polynomial of degree at most 2 A' - 1. This, together with the uniform boundedness of and 

(4.8), (4.10) and (4.11), implies (4.7) and the uniform boundedness of the quantity under the limit sign with 
respect to t. 
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We thus only have to deal with the left hand side of (4.6). We integrate (4.6) in t, integrate by parts for 
the first, term, and use Lemma 4.2 to obtain: 


-a: 


Un(x, t) 


(4.12) 


f 1 f n ^i’N-l(x,0) J 
- J ^ujv(a:,0)— y===-dx = o(l). 


Again, the difficulty is that we have only assumed the uniform boundedness of £), hence of f{uiv(xj , t )), 

not the uniform boundedness of either u N (x,t) or Jjv/(ttjv(:r,f)). We again get around this by switching 
between integrals and quadrature summations: 


/: 


u N (x,t) 




+ J)v/(ujv0 


(4.13) = jhwjy/T^ (u N ( Xj ,t) 

because we can easily verify that the integrand 


+ }{u N (x 


"*» (W 


t) 


dx 



u N (x,t) 


f ^ 7 V-l(x,t) ^ 

V %/T - x 2 ) 


+ J N f{u N {x,t)) 


{X,t)\ \ 

V vT^ U 


is a polynomial of degree at most 2N - 1. Recalling the definition of the function v Ajr (x, t) in (1.2) and that 
of <p Ax (x,t ) etc. in (2.10), we can use (4.13) to rewrite (4.12) as 




Ax 






dx dt 


(4.14) 

where 


- I' 0(*)i>a,(M) dx = 0(1) 


> Wjy/T— X 2 

(i{x) = - , X j < X < Xj+i. 

x j+l — x j 

Clearly, 0(x) is uniformly bounded and converges to 1 as Ax -> 0. By assumption, v& X {x,t) converges to 
u{x, t ) boundedly a.e. Also, (4.4) guarantees the uniform convergence of the ipN-\ related terms to the right 
limits. Using the dominant convergence theorem, taking the limit as Ax 0 in (4.14), we obtain (1.3). 
This proves that u(x,t) is a weak solution of (1.1), i.e. we have proved the following 


Proposition 4.3. If the function v& X (x,t) defined in (1.2), obtained from the solution of the Chebyshev 
collocation scheme (4.1) at the Chebyshev Gauss-Lobatto points Xj, converges almost everywhere to a 
function u(x, £), then u(x,£) is a weak solution to (1.1). 


5. Multi-Domain Legendre Methods. In this section we will discuss stable and conservative inter- 
face boundary conditions for the multi-domain Legendre method applied to equation (1.1). We assume that 
the domain -1 < x < 1 is divided into tw T o domains, and for the sake of simplicity we assume that the 
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interface point is x = 0. We will denote by u N (x, t) the numerical approximation in -1 < x < 0 and by 
v N (x,t) the solution at 0 < x < 1. The multi-domain Legendre method is given by 


(5.1) 


(5.2) 


+ = B ( U N{~l,t)) + nQ T (x) [/ + («A'(0, f)) - f + {v N (0,t))] 

+t 2 Q i (x) [f~ (u N {0,t)) - / _ (r>A-(0,t))] + SV(un), 
n) = TtQ'V) [f + (vN(0,t)) - f + (u N (0,t))] 

+t 4 Q ti {x) [f~(v N (0,t)) - /“(u^(0,0)] 
+SV(v N )+B(v N (l,t)). 


Equation (5.1) holds in the interval — 1 < x < 0, and (5.2) holds in 0 < x < 1. interpolates /(ujv) 

at the zeroes of the polynomial xQ 1 and ifj J(vn) interpolates /(uyv) at the zeroes 77 j of the polynomial 
xQ n , where 


The spectral viscosities SV(un) and SV(vjv) are of the form 


(5.3) 

(5.4) 


SV(u N ) = e 


SV(v N ) = c 


(~ 1)» +1 

AT 2 -’- 1 

(~l ) s+1 


N 2s-1 


d , . d 


S 


' d_ 

dx 


c) 

1(1 _x) s 


V N . 


At this point we stress that the results of this section are valid only for this form of spectral viscosity 
and not for the others discussed in Section 3. The reason for that will be evident in the proof. 

Finally the boundary operators B at the ends of the interval — 1 < x < 1 are left unspecified for now. 

We will also denote the scalar product (p, <?)jv = ^£=0 if p{ x ) an d Qi x ) are defined in 

[—1,0] and (p,q)N = Y,f^oP T ( r )j)<l( r )j) UJ j if P( x ) and Q( x ) are defined in [0,1], and u)j are the weights in 
the Gauss-Lobatto Legendre quadrature formula. Note that if pq is a polynomial of degree at most 2 N — 1 
defined in [—1,0] then 


( p,q)N = (p, q) 


,0 

=ij 


(x)q(x)dx. 


A similar formula holds in the interval [0, 1]. 

Our aim in this section is to show that the choice of the parameters r*, i = 1,4 that leads to linear 
stability is sufficient for proving conservation, i.e. if the numerical solution ujv(x, i), v N (x,t) converges 
boundedly a.e to functions u(x, t), v(x ) t), then the solution w defined by w(x^t) = u(x,t) if — 1 < x < 0 and 
w(x, t) = v(x, t) if 0 < x < 1 converges to the weak solution of (1.1). 

We will discuss first the stability of (5.1)-(5.2). We state 
Proposition 5.1. The boundary operators are dissipative, i.e. 

(5.5) («A r ,%iv(-M)))^ + -uJf(-l 7 t)Au n {-l>t) < 0, 

(5.6) {v N ,B(v N {l,t))) N - it>A(-M).4r’/v(M) < 0. 
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Proposition 5.1 implies that the boundary treatment at the end-points of the interval is stable. Example 
of such operators are given in [4], 

We are ready to state the stability theorem for the linear constant coefficient case. In this case there 
is no need for the spectral viscosity terms and we will ignore them. We assume that / = An where .4 is 
symmetric in equation (1.1) and in the same way / + = f~ = A~u where the eigenvalues of ,4 + are 
nonnegative and those of A~ are nonpositive. 

Theorem 5.2. Let un, vn be the solutions of (5.1)-(5.2). Define 


E(t.) = (u N (x,t),u N (x,t)) N + (v N (x,t),v N {x,t)) N , 


then 


E(t) < E( 0) 


provided that 


(5.7) 


Tl ~ 2w 0 ’ 


r 2 > 


2u)q 


7*3 < 


1 


2uj 0 

1 

n -r 3 = — , 


r 4 > , 

ZiJQ 


T-2 


T\ = . 

u/q 


Proof: The proof follows from multiplying (5.1), (5.2) by ujr,vjj and taking the scalar product. We use 
Proposition 5.1 and the following notation 


u 0 - u N ( 0, t), v 0 = v N (0 , t), oj = U 0 A ± 7l 0 , ft* ~ v 0 A ± v 0 , Jo = UqA ± Vq 


to get 




i )Q ° + ■ (n + Tsh ° + + 


(5.8) 


+ to ~ 2^ Q ° _ to + toio + (n + d-)& 


The conditions stated above for the rC s guarantee that the right hand side of (5.8) is nonpositive and 
the proof is completed. 


Remark: The Discontinuous Galerkin method applied to this problem leads to the upwinding choice T\ = 
T4 = 0, r 2 = 73 = — -j-. Another attractive choice that involves no splitting of the fluxes is T\ = r 2 = — 73 = 


We turn now to the main purpose of this section, namely the proof of convergence in the nonlinear case 
of the numerical solution to the correct entropy solution. 

We first show that the spectral superviscosity terms do not create any problems: consider a compactly 
supported (in [-1,1]) test function $(x, t) in C m [- 1, 1], m is to be specified later. 

Lemma 5.3. Let and be the Legendre interpolation polynomials of ^(x, t) in the 

intervals [-1,0] (with collocation points £j), and [0, 1] (with collocation points 7^), respectively. Then 

Jim (0jv-i,Sr(ujv)) = 0, 
lira (ip N -i,SV[v N )) = 0. 


(5.9) 

(5.10) 
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where the spectral superviscosities are defined in (5.3) and (5.4). 


Proof: Since u ^ is a polynomial of degree N it can be represented as 


u N {x) = y «*-,Ar r -P*:(2a- + 1). 


Therefore from (5.3) 


(_n«+i o cH s 

SV M = e^^ r - X { X + 1)- u N 


= -CT^rr E k ^ k + l) S ^,yvP* (2:r + 1). 


Also the test function 4>n~\ can be represented as 


&N - 1 = + 1)* 


^From the orthogonality of the Legendre polynomials it follows that 


|(0JV-i,-SV(«iv))| l)*(/ > *(2ar + 1), /\.(2a: + 1)). 

k~ 0 

We can choose m large enough such that 

C 

\$k,NUk t N\ < p 

and since and (Pk{ 2x + l),P*(2a: + 1)) are bounded, the proof is established. The proof for (5.10) is 
similar. 

■ 

It is self evident that the form of the spectral viscosity SV is crucial. In fact the factor 1 — x 2 is necessary 
in the proof. Note that 

(5.11) (4>n-i > SV(un)) = (SV(<f>N-\), wat) = (51 (<£jv_i),ujv)jv 

We basically proved that the first argument in the scalar product in the right hand side of (5.11) tends to 
zero whereas the second argument is bounded. The relation (5.11) is not true for other forms of the spectral 
viscosity where the factor 1 — x 2 does not appear. 

Lemma 5.4. Let r* satisfy (5.8), then 

,, 0u N d(f > N - 1 , ,, dv N f . d$N-\. 

i> ~of) N “ (/(«iv)» )n + Wv-i, ” (/Wv)> )at 

(5.12) = (<^v_i,5I (u^)) + (tl>N-hSV(vN)) 


Proof: Taking the scalar product of equation (5.1) with (f>N-\ and (5.2) with Vw-i and denoting by /( 0, t) = 
/o, for all the quantities one gets 

{<f>N- 1, ”^~)A' T (<t>N-U ^ — r 1^0 (/ + ( w o) — />o)) ^0 
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(5.13) 


+7200 (/ («o) - / (»’o)) W 0 
+(0jv- i, 5 V'(wjv))jv, 

W>JV-1, + (0AT-1, = T 30O (/ + (7’o) - / + («0 )) W 0 


(5.14) 

Here we used the fact that 0(— 1, t) = 0(1, £) — 0. 
Now, 


7400 (/ (i’o) - / («o)) w 0 
+ (0 N — 1 , 5 V (t’n ) ) JV ) • 


(0JV-1, 


dljjM 

dx 


(5.15) 

and by the same token 

(5.16) 


\ _ <i oiy(u N ) 

)N “ (0AT-1, ) 

= 0o/(mo) - (/^/(«jv), 1 ) 

= 4>of(Uo)-(fM,^ 1 )N 


dljJf(v N ), 


dv>N - 1 


(0yv-i , - - j; — )n = -MM - (f{v N ), dx 

Using (5.15)-(5.16) in (5.13)-(5.14) one gets 
' du j\ T 


)n- 


<t>N- 


dt 


(/(«*). + 


( s dvN \ 

{ VN - U ~df) ~ 


/ rf \ dtyN—\ \ 

{fiv A r ), dx )n 


- (/ + (flo) “ f + M) [t\Vo - T3UJ0 - 1 ] <£o + (/ (w 0 ) - / (uo)) [T2UQ - t 4 u ? 0 - 1] <f>0 
+ {(f> N _i,SV{u N )) + {iI>n-uSV(vn)) . 


Taking (5.8) into account, the lemma is proven. 


We integrate now (5.12) with respect to time to get 



= “ (Uw{t = O),0AT-l(t = 0)) - {v N (t = O),'0N-l(* = 0)) . 


We now use Lemma 3.3 in Section 3 to convert ujy, vj\< to u& x which are defined in (1.2) as the 
piecewise polynomials having the values of ujv, at the grid points. Combining then with Lemma 5.3 it 
follows: 

Theorem 5.5. Let u/v and be the multi-domain Legendre approximation (5.1)-(5.2) to (1.1). Assume 
that the functions u& x and v& x defined in (1.2) converge boundedly a.e., then the limit function is a weak 
solution of (1.1). 
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